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ABSTRACT 

We present radial velocity measurements of two stars observed as part of the Lick Subgiants Planet Search 
and the Keck N2K survey. Variations in the radial velocities of both stars reveal the presence of Jupiter-mass 
exoplanets in highly eccentric orbits. HD 16175 is a GO subgiant from the Lick Subgiants Planet Search, orbited 
by a planet having a minimum mass of 4.4 Mj up , in an eccentric (e = 0.59), 2.71 yr orbit. HD 96167 is a G5 
subgiant from the N2K ("Next 2000") program at Keck Observatory, orbited by a planet having a minimum 
mass of 0.68 Mj up , in an eccentric (e = 0.71), 1.366 yr orbit. Both stars are relatively massive (M* = 1.3 M e ) 
and are very metal rich ([Fe/H] > +0.3). We describe our methods for measuring the stars' radial velocity 
variations and photometric stability. 

Subject headings: extrasolar planets — stars: individual (HD 16175, HD 96167) — stars: planetary systems — 
techniques: radial velocities 



1. INTRODUCTION 

Over 300 exoplanets have been detected to date, 9 a quan- 
tity that is making it possible to characterize statistically sig- 
nificant trends that shed light on planet formation processes. 
Trends in planet occurrence as a function of host star mass 
and metallicity have emerged, and the distribution of orbital 
eccentricities is still being un covered. For a r eview of ex- 
oplane t detections , see, e.g., Ijohnsonl d2009l) : Marc y et al.l 
d2008h : ludrvet al.ld2007b . 

Both the planet occurrence rate and the total planetary 
mass in a system are proving to be positively correlated 



with t he host star's mass (Johns on et al . 2007; lLovis & M ayor 
2007b . Our understanding of that planet occurrence rate at 
masses of M+ > 1 .3 M Q comes from surveys of evolved stars. 
Intermediate-mass stars cool and slow their rotation as they 
evolve off the main sequence, thereby developing the copi- 
ous narrow metal lines that facilitate precision radial veloc- 
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ity monitoring. Since giant stars tend to have considerable 
photospheric jitter (introducing errors of ~20ms _1 ), sub- 
giants, which have jitters of ~ 5 ms , represent an evolution- 
ary sweet spot for detecting planets around intermediate-mass 
stars. 

In addition to correlating with stellar mass, the planet oc- 
currence rate is observed to correlate with stellar metallicit y 
dFischer & Valentill2005l ISantos et al.ll2004l lGonzalez|[l997l) . 



The results of Doppler surveys are corroborated by theoreti- 
cal predictions that more metals in a protoplanetary di sk leads 
to more prolific planet formation dlda & Linl 12004 ). This 
trend informs target selection for s ome exoplanet searc hes, 
e.g., the N2K program ("Next 2000:' lFischer et aljj2005|) _and 



the E LODIE Metallicity-Biased Planet Search dda Silva et al 
2006), which specifically target metal-rich stars. Such 



searches allow planet hunters to maximize exoplanet discov- 
eries given finite telescope resources. 

Exoplanet systems exhibit a broad range of eccentricities. 
For the 267 well-characterized planets within 200 pc (most of 
the known exoplanets), those in non-tidally-circularized orbits 
(P > 5 days) span a range of eccentricities from 0.00 to 0.93, 
with a median eccentricity of 0.23. The broad range of ec- 
centricities is surprising, since it was previously expected that 
such objects would have orbits resembling the circular orbits 
of the solar system gas giants. The origin of exoplanet eccen- 
tricities is one of the major outstanding problems in the study 
of exoplanets. The discovery of additional planetary systems, 
especially those with large eccentricity, will bring the shape 
of the eccentricity distribution into better focus and guide the- 
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We present here two new Jupiter-mass exoplanets detected 
with the radial velocity method. Both stars are metal-rich 
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([Fe/H] > +0.3), intermediate-mass (Af* = 1.3 M Q ), G-type 
subgiants, and therefore contribute to our understanding of 
the planet occurrence rate as a function of mass and metal- 
licity. HP 16175 is p art of the Lick Subgiants Planet Search 
(I Johnson et all 120061) and HP 96167 is part of the N2K pro- 
gram at Keck Observatory (iFischer et al.ll2005l) . Because the 
planets presented here are from different programs at differ- 
ent observatories, we describe the observational, data analy- 
sis, and Keplerian techniques for each star in its own section. 
HP 16175 is presented in §E]and HP 96167 in §[3] We review 
the photometric stability of the stars in §|4]and summarize our 
results in §|5] 

2. HD 16175 

2.1. Observations and Radial Velocity Measurements 

We began Poppler monitoring of HP 16175 (HIP 12191, 
BP+41 496) at Lick Observatory in 2004 November as part 
of the Lick Subgiants Planet Search, a program comprising 
159 evolved intermediate-mass stars, chosen to have color 
index 0.55 < B — V < 1.00 mag, V-magnitude brighter than 
7.6 mag, and declination d > —20°. In order to choose 
evolved stars, we require an absolute visual magnitude My 
more than 1 .0 mag above the main se quence. For th is purpose, 
we use the main sequence defined in Wright| [2005. which is a 
ninth-order polynomial fit to the Hipparcos dwarf stars within 
60 pc. In order to ensure selection of subgiants rather than gi- 
ants, we required My fainter than 1.0 mag. The color limits 
select stars for which radial velocity monitoring is feasible, 
while the V-magnitude and declination limits allow the sam- 
ple to be observed with the 3 m Shane telescope and with the 
0.6 m Coude Auxiliary Telescope (CAT), both of which feed 
the L ick O bservatory Hamilton echelle spectrograph dVogt 
19871) . See Ijohnson etail J2OO6) for a complete description 
of the Lick Subgiants Planet Search and its selection criteria. 

We use the Hamilton spectrograph at a setting that provides 
resolution R w 50000 at X = 5500 A. On the CAT, typical ex- 
posure times were 60 minutes, yielding signal to noise ratio 
(SNR) of ~ 50 at 5500 A. On the Shane, shorter 10-minute ex- 
posures generally yielded higher SNR ~ 175. The lion's share 
of the data presented here were collected with the CAT. The 
procedure for ob taining radia l veloc ities from Hamilton spec- 
tra is outlined in iButler et alJ (fl996). In short, a temperature- 
controlled Pyrex cell, filled with gaseous molecular iodine, 
is placed before the spectrometer's entrance slit. The iodine 
imprints a dense set of absorption lines between X ~ 5000 A 
and X ~ 6400 A, which are calibrated to provide an accurate 
wavelength scale for each observation, and which also serve 
to provide information about the spectrometer's instrumental 
response dMarcv & Bufleril 1992b . 

The Poppler measurement method described above tradi- 
tionally includes the use of a template spectrum, an iodine- 
free, high-S/N observation relative to which the Poppler 
radial velocity measurements are made. In the case of 
HP 16175, we at first measured the Poppler shifts relative to 
a synthetic "morphed" template spectrum. The procedure for 



TABLE 1 
Radial Velocities for HD 16175. 
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Uncertainty 
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4.68 


14425.847 
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14426.798 
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14525.640 
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14548.644 


26.06 
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14675.985 
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-22.00 


6.32 


14699.954 


-37.38 


11.05 


14734.937 
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5.39 


14748.784 


-82.01 


6.90 


14767.935 


-94.18 


5.89 


14846.680 


38.50 


5.49 


14858.657 


43.15 


7.74 


14904.639 


78.14 


7.19 



creatin g the morphed template is described in Ijohnson et al 
(2006). Once the star showed the potential for a planetary 
companion, we obtained a traditional template observation. 
The observed template is also used for the LTE analysis de- 
scribed in § 12.21 

The radial velocities, measured in approximately 700 2 A 
segments of each echelle spectrum, are averaged to determine 
the radial velocity of a star in a given observation. Table Q] 
contains the velocities and their associated uncertainties for 
HP 16175, from 2004 November to 2009 March. We esti- 
mate those uncertainties to be the standard deviation from the 
mean of the 700 individual segment velocities. Typical inter- 
nal uncertainties for HP 16175 are ^6 ms -1 . The internal 
uncertainties in Table [T] are added in quadrature to a 5 ms 
error due to stellar j itter (a typical figure for s ubgiant stars; see 
Fischer et alj|2003l and lJohnson et al.ll2.007b . and those com- 
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bined values are used in the least-squares Keplerian analysis 
described in § 12.31 



TABLE 2 
Stellar Parameters. 



2.2. Stellar Properties 

HD 16175 is a GO subgiant in the foreground of galactic 
open cluster M34 (NGC 1039), which is 440 pc farther away. 
In our analysis we found several different literature V-band 
magnitudes for HD 16175, so we adopted for ou r analysis here 
an ave rage of the values from three sources: the iTere shchenkol 
(2001) spectrop hotometric standards list, the Hipparcos cata- 



log dESAlll997l) . and the Sky2000 catalog. 10 Our adopted 



values were V = 7.29 ± 0.03 and B — V = 0.63 ± 0.02. 

The Hipparcos parallax distance for HD 16175 is 60 ±3 pc, 
which, when taken together with V from the preceding para- 
graph, implies an absolute visual magnitude My = 3.40. Its 
My places HD 16 175 1.3 mag above the Hipparcos main se- 
quence defined in lWrightl (2005, a value that fulfills the selec- 
tion criterion described in § 12. 11 1. HD 16175's high metallicity 
(measured in the next paragraph) may account for part of its 
distance above the main sequence, but some of that 1.3 mag 
difference is likely due to its evolution, even though it may 
not yet have exhausted its core hydrogen. 

HD 16175 is metal rich, with [Fe/H] = +0.39 ±0.06, and 
has an effective temperature T e s — 6080 ± 70 K. Metallicity 
and temperature come from our LTE spec tral synthesis analy- 
sis wi th Spectroscopy Made Easy (SME; IValenti & Piskunov 
1996b. In general we follow the procedure outlined in 
Valenti & Fischer] ([2005 ). The LTE analysis also yielded 
surface gravity logg = 4.4 ± 0.2 and projected rotational 
velocity vsin/ = 4.8 ± 0.5 kms~'. The uncertainties of 
the SME-derived stellar properties ([Fe/H], T^, logg, and 
v sin Q were determin ed f ollowing the procedure de scribed in 
Johnson et all d2008l) and I Valenti & Fischerl d2005l) . The full 
set of stellar properties and their uncertainties (in parentheses) 
appear in Table [2] 

The stellar luminosity L± = 3.49 ± 0.12 Lp, was calcu- 



lated using the relation in IValenti & Fischerl (120051). based 



on My and a bolometric correction from VandenB erg & Clem 
(2003). We follow the same procedure here, adopting a bolo- 
metric correction of —0.035 mag. The stellar radius — 
1.68 ± 0.04 R was calculated with the Stefan-Boltzmann 
law, using L* given and T e s from the LTE analysis. We 
then calculated the stellar mass M* = 1 . 29 ±0 .09 M Q and 



age 4.0 ± 1 .0 Gyr using the iGirardi et al.l 020021) stellar inte- 
rior models, based on the Hipparcos-deiived My, the B — V 
color, and the LTE-derived metallicity. The uncertainties for 
these derived values are propagated from the uncertainties on 
the input values. 

The chromospheric activity parameters Shk and logR ! HK , 
based on the Ca II H and K emission features, are mea- 



sured according to the procedure in Wri ght et al.l (120041) . The 
full listing of stellar characteristics for HD 16175 appears in 
Table |2] 



HD96167 a 



8.09 
3.46 
0.73 
84(9) 
+0.34 (0.06) 
5770 (70) 
3.8 (0.5) 
4.0 (0.2) 
1.31 (0.09) 
1.86 (0.07) 
3.4 (0.2) 
-0.081 
3.8 (1.0) 
0.14 
-5.16 



Parameter HD16175 a 

V (mag) 7.29 

My (mag) 3.40 

B-V(mag) 0.63 

Distance (pc) 60(3) 

[Fe/H] +0.39 (0.06) 

r eff (K) 6080(70) 

vsin((kms-') 4.8(0.5) 

logg 4.4 (0.2) 

W. (M ) 1.29 (0.09) 

«. (R. 1 1.68 (0.04) 

/.. (I.. ) 3.49 (0.12) 

Bolometric Corr. (mag) —0.035 

Age (Gyr) 4.0(1.0) 

■Shk 0.17 

logR^K 496 

1 Uncertainty given in parentheses. 



2.3. Keplerian Fit 

We search for the best fitting Keplerian orbital solution to 
the radial velocity time series with a Levenberg-Marquardt 
least-squares minimization. Beginning with a periodogram 
to determine an initial guess for the orbital parameters, our 
Keplerian fitter searches a grid of initial period guesses to 
determine the orbital solution. We estimate the uncertain- 
ties in those orbital parameters with a bootstrap Monte Carlo 
method: we subtract the best-fit Keplerian from the measured 
velocities, scramble the residuals and add them back to the 
original measurements, then obtain a new set of orbital pa- 
rameters. We repeat this process for 1000 trials and adopt the 
standard deviations of the parameters from all trials as the for- 
mal uncertainties. For a more in-depth description of the Ke- 
pleri an fitting and uncer tainty determination procedures, see, 



Marcv et al. (2005). 



e.g., 

For HD 16175, the best fitting Keplerian orbital solution 
has a period of P = 990 ± 20 days (or 2.71 ± 0.05 yr). It also 
has eccentricity e = 0.59 ± 0. 1 1 and velocity semi-amplitude 
K = 94±11 mr 1 , The rms scatter of the data about the fit 
is 9.2 ms" 1 , and the reduced ^/Xv i s 116. Based on our 
stellar mass estimate in § 12.21 of M* = 1.29 M Q , the planet 
HD 16175b has a semimajor axis a = 2.1 AU and a minimum 
planet mass Mp sin i = 4.4 Mj up . The full list of orbital param- 
eters appears in Table[3] with their associated uncertainties in 
parentheses. 

The velocities from Table [T] are plotted in Figure [T] The 
error bars on the plot represent the quadrature addition of the 
internal measurement errors from Tableland 5 m s~ 1 of stel- 
lar jitter (see § 12. U . The best-fit orbital solution is plotted 
with a dashed line and the orbital parameters are printed in 
the plot. The photometric stability and transit probability for 
HD 16175 b are discussed in §0] 



1 Vizier Online Data Catalog, 5109 (Myers et al. 2001) 
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TABLE 3 
Orbital Parameters. 



Parameter HD 16175 b a HD 96 1 67 b 

P (d) 990 (20) 

/.,'• (JD) 2452820 (30) 

e 0.59(0.11) 

K(ms , 94(11) 

Cfl(deg) 222 (9) 

Mpsini (M Jup ) 4.4(0.34) 

a(AU) 2.1 (0.08) 

Fit RMS (ms" 1 ) 9.2 

Jitter (ms~') 5.0 

Reduced- \f$ 1.16 

AU S 44 



498.9(1.0) 
2453057 (5) 
0.71 (0.04) 
20.8 (1.5) 
285 (7) 
0.68 (0.18) 
1.3 (0.07) 
4.6 
5.0 
0.88 
47 



Uncertainty given in parentheses. Time of periastron pas- 



3. HD 96167 

3.1. Observations and Radial Velocity Measurements 

Doppler monitoring of HD 96167 (HIP 54195, 
BD -09 3201) began at Keck Observatory in 2004 Jan- 
uary as part of the N2K consortium observing program. 
N2K targeted 2000 stars within HOpc with the goal of 
detecting high transit probability hot Jupiters. The stars were 
chosen to have 0.4 < B — V < 1.2 mag and V brighter than 
10.5 mag. The targets were also chosen to be metal rich, 
with [Fe/H] of +0-1 o r higher, as determined photometrically 
(Ammonsetal. 2006) and with low-resolution spectroscopy 
Robinson et al . 2006). Observations each star were collected 



in a high-cadence time series intended to detect short-period 
planets. F or a full descr i ption of the N2K consortium and its 
goals, see Fischer et al.l (120051) . The main phase of the N2K 



program has finished, but several stars that showed significant 
velocity scatter have been monitored further. A number of 
longer period planets have been detected around these stars, 
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FIG. 1. — Radial velocity time series for HD 16175, measured at Lick 
Observatory. The error bars represent a quadrature sum of the internal mea- 
surement error (listed in TablefTJ and 5 ms~' stellar jitter (see j) |2.3| for more 
detail). The dashed line represents the best Keplerian orbital fit to the data 
for HD 16175 b. The parameters of that fit are printed in the plot inset. The 
residuals between the Keplerian fit and the data appear in the lower panel of 
the plot. 



including the planet orbiting HD 96167 announced here. 

HD 9 6167 was observ ed with the HIRES echelle spectro- 
graph (lVogt et al.lll994l) on Keck I. The spectra have R rj 
70000 at X = 5500 A. Typical exposure times were 2 minutes, 
yielding SNR ~ 250 at 5000 A. The iodine cell procedure for 
measuring radial velocities in Keck spectra is the same as for 
Lick spectra, described in § 12.11 As with the Lick Subgiants 
Planet Search, we initially observed the N2K stars relative to a 
morphed synthetic template spectrum (see § 12. U . taking a tra- 
ditional template exposure once a planet candidate emerged. 

The full list of measured velocities and their associated 
internal uncertainties (estimated as described in § 12. lb ap- 
pears as Table [4] The observations span the period from 
2004 January to 2009 January. Typical unc ertainties are 
l , and a jitter un certainty of 5 ms _1 dFischer et al 



1 ms 



2003; John son et al.ll2007l) is added in quadrature before the 
least-squares Keplerian fit is carried out (§ 13. 31 . 

3.2. Stellar Properties 

HD 96167 is a G5 subgiant star. In contrast to HD 16175, 
literature V-band magnitudes for HD 96167 were consistent 
with Hipparcos, so we adopt the Hipparcos values here: V — 
8.09 andB-V = 0.73. 

HD 96167 has a Hipparcos parallax-based distance of 84 ± 
9 pc. The parallax distance and V-magnitude, when taken 
together, yield My = 3.46, plac ing it 1.5 ma g above the 
main sequence as defined by the IWrighll (120051) polynomial. 
HD96167's high metallicity (measured in next paragraph) 
may account for part of its distance above the main sequence, 
but its My places it far enough along the evolutionary track 
for 1.3-M stars that it has likely exhausted its core hydrogen 
and evolved into a subgiant. 

Fitting with its inclusion in the N2K stellar sample, 
HD 96167 is metal rich, with [Fe/H] = +0.34 ±0.06, and 
has a near-solar temperature T e g = 5770 ± 70 K. It has 
logg = 4.0±0.2 and vsin; = 3.8±0.5 kms -1 . Its mass is 
1.31 ±0.09 M and its age is 3.8 ± 1.0 Gyr; it has a lumi- 
nosity L* of 3.4 ±0.2 L and radius of 1.86 ± 0.07 R . 
Its Ca II H and K chromospheric activity is relatively quiet, 
with S H k = 0.14 and \ogR^ K = -5-16. As with HD 16175, 
the [Fe/H], T e f[, logg, and vsin; come from an LTE spectral 
synthesis using SME. The procedures for determining L*, 
M*, age, Shk, and logR' HK and their associated uncertainties 
are identical to those described in § 12.21 The full set of stellar 
parameters is listed in Table[2] 

3.3. Keplerian Fit 

The procedure for determining the best fitting Keplerian or- 
bital solution for the velocities listed in Table [4] is the same 
as described in § 12.31 The best fitting orbital solution has 
P = 498.9 ± 1 .0 days (or 1 .366 ± 0.003 yr). The eccentric- 
ity e is 0.71 ±0.04 and the velocity semi-amplitude K is 
20.8 ± 1.5 ms , The residuals from the best fit have an 
rms scatter of 4.6 ms~'; the fit has a reduced ^/Xv °f 0.88. 
Based on the LTE-determined mass estimate (§ 13. 2b of M* = 
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TABLE 4 
Radial Velocities for HD 96 167. 



Radial Velocity Uncertainty 



JD- 2,440,000 



(ms" 1 ) 



(mi" 1 ) 



13015.110 
13016.127 
13017.114 
13044.149 
13046.017 
13069.005 
13073.938 
13153.798 
13195.752 
13196.783 
13197.766 
13369.127 
13370.114 
13398.012 
13398.943 
13400.962 
13424.993 
13425.941 
13426.978 
13428.094 
13431.095 
13478.938 
13480.838 
13483.769 
13546.752 
13724.121 
13748.102 
13753.044 
13754.051 
13775.991 
13776.990 
13806.985 
14131.028 
14216.853 
14248.838 
14428.146 
14429.120 
14493.008 
14547.913 
14548.850 
14602.805 
14603.795 
14639.781 
14806.156 
14811.120 
14839.123 
14847.065 



-18.65 
-22.75 
-21.43 
-12.93 
-17.81 
19.73 
21.05 
0.74 
-9.08 
4.47 
3.06 
-8.74 
7.48 
-13.93 
-13.28 
-15.04 
-11.88 
-16.44 
-15.67 
-13.57 
-13.24 
-15.86 
-14.61 
-18.05 
-9.70 
2.03 
-5.92 
-1.56 
-7.94 
-0.60 
-4.41 
-9.75 
8.90 
4.74 
1.79 
-10.12 
-13.51 
-20.36 
-2.72 
-6.34 
18.48 
12.53 
6.98 
4.27 
-6.34 
-8.13 
-15.94 



1.99 
1.82 
1.98 
1.48 
1.62 
2.12 
2.32 
1.70 
1.51 
1.53 
1.56 
1.18 
1.33 
1.14 
1.11 
1.10 
1.03 
1.16 
1.20 
1.14 
1.10 
1.20 
1.11 
1.16 
1.62 
1.09 
1.37 
0.82 
0.94 
1.24 
1.58 
1.60 
1.50 
1.63 
1.41 
1.05 
1.15 
1.17 
1.10 
1.96 
1.40 
1.46 
1.34 
1.66 
1.25 
1.28 
1.31 



1.31 M©, the planet HD 96167 b has semi-major axis a = 1.3 
and minimum mass Mpsini = 0.68 Mj up . 

The HD 96167 radial velocities listed in Table |4] are plot- 
ted in Figure[2] As with HD 16175, the error bars represent 
the quadrature sum of the internal measurement uncertainties 
from Table 2] and stellar jitter of 5 ms -1 . The best-fit orbital 
solution appears as the dashed line in Figure |2] The pho- 
tometric stability and transit probability for HD 96 167 b are 
discussed in §|4] 

4. PHOTOMETRIC MEASUREMENTS OF HD 16175 AND HD 96167 

We have obtained three years of high-precision differen- 
tial photometry of both HD 16175 and HD 96167 with the 
T12 0.8 m automatic photometric telescope (APT) at Fairborn 
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FIG. 2. — Radial velocity time series for HD 96167, measured at Keck 
Observatory. The error bars represent a quadrature sum of the internal mea- 
surement error (listed in Tablef4) and 5 ms -1 stellar jitter (see jj |3.3l for more 
detail). The dashed line represents the best Keplerian orbital fit to the data 
for HD 96167 b. The parameters of that fit are printed in the plot inset. The 
residuals between the Keplerian fit and the data appear in the lower panel of 
the plot. 



Observatory. The APT detects short-term, low-amplitude 
brightness variability in solar-type stars due to rotational 
modulation in the visibility of photospheric starspots (e.g., 
Henry et al.lll995l) . as well as lon ger-term variations associ- 
ated with stellar magnetic cycles (Henry| |l999l) . We obtain 
photometric observations to establish whether observed radial 
velocity variations in a star are due to reflex motion caused by 
a planet ary companion or due to the effects of stellar activ- 
ity (e.g.. bueloz et alJbOQlt IPaulson et al.ll2004 . Photomet- 
ric observations can also lead to the detection of planetary 
transits and the direct determination of planetary radii, as in 
Winn et alj d2008b . 

We acquired 181 observations of HD 16175 between 2006 
December and 2009 February and 319 observations of 
HD 96167 between 2004 November and 2007 May. The T12 
APT uses two temperature-stabilized EMI 9124QB photo- 
multiplier tubes to measure photon count rates simultaneously 
through Stromgren b and y filters. On a given night, the tele- 
scope observes each target star and three nearby comparison 
stars, along with measures of the dark count rate and sky 
brightness in the vicinity of each star. Designating the com- 
parison stars as A, B, and C, and the target star as T, the ob- 
serving sequence is as follows: dark, A, B, C, T, A, skyA, B, 
sky B , C, sky c , T, sky T , A, B, C, T. For HD 16175, the com- 
parison stars A, B, and C were HD 16176, HD 14095, and 
HD 14064, respectively; for HD 96167, the comparison stars 
were HD 96220, HD 95938, and HD 94206. A diaphragm size 
of 45" and an integration time of 20 seconds were used for all 
integrations. 

The measurements in each sequence were reduced to form 
three independent measures of the six differential magnitudes 
T - A, T - B, T - C, C - A, C - B, and B - A. These differen- 
tial magnitudes were corrected for extinction and transformed 
to the standard Stromgren photometric system. To increase 
the S/N, the data from the b and y passbands were averaged 
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TABLE 5 

Summary of the Ensemble Photometric Observations 





Star 
(1) 


APT 

(2) 


Date Range 
(JD-2,440,000) 
(3) 


(4) 


Yearly Mean 
(mag) 
(5) 


a 

(mag) 
(6) 


(mag) 
(7) 


HD 16175 




T12 


14080-14172 


47 


0.41236 


0.00123 


0.00018 








14370-14533 


86 


0.41232 


0.00160 


0.00017 








14730-14881 


48 


0.41079 


0.00124 


0.00018 


HD 96167 




T12 


13330-13511 


130 


0.98318 


0.00136 


0.00012 








13690-13886 


124 


0.98380 


0.00144 


0.00013 








14058-14255 


65 


0.98287 


0.00107 


0.00013 



to create "(b+y)/2" magnitudes. After passing quality con- 
trol tests, the three independent measures of each differential 
magnitude were combined to give a single mean data point per 
complete sequence for each of the six differential magnitudes. 

To improve our measurement precision still further, we av- 
eraged the three mean T — A, T B, and T C differential 
magnitudes from each sequence into a single value represent- 
ing the difference in brightness between the program star and 
the mean of the three comparison stars T — (A + B + C)/3, 
which we refer to as the ensemble mean. This helps to av- 
erage out any subtle brightness variations in the three com- 
parison stars. The ensemble means for the three observing 
seasons of both stars are summarized in Table [5] and plotted 
as open circles in Figure[3] 

Column 5 gives the yearly mean T — (A + B + C)/3 dif- 
ferential magnitudes in (b + y)/2 for the three observing sea- 
sons. Column 6 gives the standard deviations, a, of these 
ensemble differential magnitudes from the yearly mean and 
so provide a measure of night-to-night brightness variations. 
Typical standard deviations for constant stars fall in the range 
0.0012-0.0017 mag for this telescope. All six of the stan- 
dard deviations in Table[5]fall within this range. Periodogram 
analyses of the three observing seasons for each star found 
no significant periodicity in the range of 1-100 days. Both 
HD 16175 and HD96167 are constant from night to night to 
the limit of our measurement precision; stellar activity due to 
starspots should have no significant effect on the measured 
radial velocities. 

The standard deviations of the yearly means are given in 
column 7, computed from the a in column 6 divided by the 
square root of the number of observations in column 4; these 
values provide an estimate of the precision of each yearly 
mean. The observed yearly means are plotted in Figure [3] as 
filled circles connected with solid line segments. For both 
HD 16175 and HD 96167, the observed scatter in the yearly 
means is somewhat higher than the estimated precision of 
the means, by factors of 5.0 and 3.7, respectively, imply- 
ing that there are slig ht year-to-yea r variations in the bright- 
ness of the two stars. iHenrvl ([1999) demonstrates that subtle 
brightness variability can be measured for solar-type stars and 
also demonstrates observationally that yearly means of con- 
stant stars can be m easured to 0.0002 mag with the APTs. 
Wright et al.1 ( I2008I) includes another example of this mea- 



surement precision for comparison stars. Thus, the long-term 
variation over a range of 0.0016 mag in HD 16175 is likely to 
be real. The range for HD 96167 is less than a millimagnitude 
over the three observing seasons and so is harder to establish 
as real without a longer time series. As subgiants with very 
low chromospheric activity levels (5hk and log R' HK in Table 
|2|, both stars are expected to have little brightness variation. 
Given the moderately long orbital periods and the current pre- 
cision of the orbits, a transit search is premature. 

5. SUMMARY 

We present here two Jovian-mass planets in relatively long 
period, eccentric orbits around intermediate-mass stars with 
supersolar metallicity. HD 16175 is a 1.29-M GO subgiant 
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FIG. 3. — Nightly mean differential magnitudes of both stars (open circles). 
The yearly means are plotted as filled circles connected with straight line 
segments. The uncertainties of the yearly means, as given in column 7 of 
Table [5] approximately match the size of the plotted symbols. The dashed 
line in each panel represents the mean of the three yearly means. The total 
range in magnitudes of the three yearly means is given in the lower left corner 
of each panel. The standard deviation of the yearly means from the mean of 
the means (dashed lines) is given in the lower right corner. 
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from the Lick Subgiants Planet Search at Lick Observatory, 
orbited by a planet with a minimum mass of 4.4 M Jup in a 2.71 
yr, eccentric (e = 0.59) orbit. HD 96167 is a 1 .31-M© G5 sub- 
giant from the N2K program at Keck Observatory, orbited by 
a planet with a minimum mass of 0.68 Mj up in a 1.366 yr, ec- 
centric (e = 0.71) orbit. The relatively high metallicities and 
masses of the stars in question make these planets particularly 
interesting additions to the exoplanet menagerie. 

Approximately half of the planets orbiting subgiants de- 
tected to date are in orbits longer than a year. HD 16175 b 
and HD 96167 b are additions to that category, with peri- 
ods of 2.71 ± 0.05 and 1 .366 ± 0.003 years, respectively. Of 
the 267 known Doppler exoplanets with well constrained 
orbits within 200 pc, 23 are in orbits more eccentric than 
HD 16175 b (e = 0.59) and only 12 are in orbits more eccen- 
tric than HD 96167b (e = 0.71). Only eight of the known 
Doppler planets orbit stars more metal rich than HD 16175 
([Fe/H]= +0.39), while only 23 orbit stars more metal rich 
than HD 96167 ([Fe/H]= +0.34). The planets fall into the up- 
per reaches of three emerging trends in exoplanet distribution 
(eccentricity, host star mass, and metallicity), thereby con- 
tributing to our understanding of the nature of exoplanets and 
of planet formation. 

We are indebted to the many observers who collected spec- 
tra of these stars. For HD 16175, we thank the observers 



Kelsey Clubb, Julia Kregenow, Joshua Peek, Karin Sand- 
strom, and Julien Spronck. For HD 96167, we thank the 
observers Gaspar Bakos, R. Paul Butler, Chris McCarthy, 
Guillermo Torres, Steven S. Vogt, and Joshua Winn. We 
gratefully acknowledge the dedication and support of the Lick 
and Keck Observatory staffs, in particular Tony Misch for 
support with the Hamilton Spectrograph and Grant Hill and 
Scott Dahm for support with HIRES. JAJ is an NSF Astron- 
omy and Astrophysics Postdoctoral Fellow with support from 
the NSF grant AST-0702821. DAF acknowledges research 
support from NASA grant NNX08AF42G. GWH acknowl- 
edges that Automated Astronomy at Tennessee State Univer- 
sity has been supported by NASA and NSF as well as Ten- 
nessee State University and the State of Tennessee through 
its Centers of Excellence program. We thank the NASA Ex- 
oplanet Science Institute (NExScI) for support through the 
KPDA program. We thank the NASA, NOAO, and UCO/Lick 
telescope assignment committees for allocations of telescope 
time. The authors extend thanks to those of Hawaiian ances- 
try on whose sacred mountain of Mauna Kea we are privileged 
to be guests. Without their kind hospitality, the Keck obser- 
vations presented here would not have been possible. This 
research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France, and of NASA's Astrophysics Data 
System Bibliographic Services. 



REFERENCES 



Ammons, S. M., Robinson, S. E., Strader, J., Laughlin, G., Fischer, D., & 

Wolf, A. 2006, ApJ, 638, 1004 
Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C, Dosanjh, P., & Vogt, 

S. S. 1996, PASP, 108, 500 
Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 

580 

da Silva, R., Udry, S., Bouchy, F., Mayor, M., Moutou, C, Pont, F, Queloz, 

D., Santos, N. C, Segransan, D., & Zucker, S. 2006, A&A, 446, 717 
ESA. 1997, The Hipparcos and Tycho catalogues, Garching 
Fischer, D. A., Laughlin, G, Butler, P., Marcy, G., Johnson, J., Henry, G., 
Valenti, J., Vogt, S., Ammons, M., Robinson, S., Spear, G., Strader, J., 
Driscoll, P., Fuller, A., Johnson, T., Manrao, E., McCarthy, C, Mufioz, 
M., Tah, K. L., Wright, J., Ida, S., Sato, B., Toyota, E., & Minniti, D. 2005, 
ApJ, 620, 481 

Fischer, D. A., Marcy, G. W., Butler, R. P., Vogt, S. S., Henry, G. W., 
Pourbaix, D., Walp, B., Misch, A. A., & Wright, J. T. 2003, ApJ, 586, 
1394 

Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 
Ford, E. B., & Rasio, F. A. 2008, ApJ, 686, 621 

Girardi, L., Bertelli, G., Bressan, A., Chiosi, C, Groenewegen, M. A. T., 

Marigo, P., Salasnich, B., & Weiss, A. 2002, A&A, 391, 195 
Gonzalez, G. 1997, MNRAS, 285, 403 
Henry, G. W. 1999, PASP, 111, 845 

Henry, G. W., Fekel, F. C, & Hall, D. S. 1995, AJ, 1 10, 2926 
Ida, S., & Lin, D. N. C. 2004, ApJ, 616, 567 
Johnson, J. A. 2009, ArXiv e-prints 

Johnson, J. A., Butler, R. P., Marcy, G. W., Fischer, D. A., Vogt, S. S., Wright, 

J. T, & Peek, K. M. G. 2007, ApJ, 670, 833 
Johnson, J. A., Marcy, G. W., Fischer, D. A., Henry, G. W., Wright, J. T, 

Isaacson, H, & McCarthy, C. 2006, ApJ, 652, 1724 
Johnson, J. A., Marcy, G. W., Fischer, D. A., Wright, J. T., Reffert, S., 

Kregenow, J. M., Williams, P. K. G., & Peek, K. M. G. 2008, ApJ, 675, 

784 

Juric, M., & Tremaine, S. 2008, ApJ, 686, 603 
Lovis, C, & Mayor, M. 2007, A&A, 472, 657 
Marcy, G. W., & Butler, R. P. 1992, PASP, 104, 270 

Marcy, G. W., Butler, R. P., Vogt, S. S., Fischer, D. A., Henry, G. W., 
Laughlin, G., Wright, J. T, & Johnson, J. A. 2005, ApJ, 619, 570 



Marcy, G. W., Butler, R. P., Vogt, S. S., Fischer, D. A., Wright, J. T, Johnson, 
J. A., Tinney, C. G., Jones, H. R. A., Carter, B. D., Bailey, J., O'Toole, S. J., 
& Upadhyay, S. 2008, Physica Scripta Volume T, 130, 014001 

Myers, J. R., Sande, C. B., Miller, A. C, Warren, Jr., W. H, & Tracewell, 
D. A. 2001, VizieR Online Data Catalog, 5109, 

Paulson, D. B., Saar, S. H, Cochran, W. D., & Henry, G. W. 2004, AJ, 127, 
1644 

Queloz, D., Henry, G. W., Sivan, J. P., Baliunas, S. L., Beuzit, J. L., Donahue, 
R. A., Mayor, M., Naef, D., Perrier, C, & Udry, S. 2001, A&A, 379, 279 

Robinson, S. E., Strader, J., Ammons, S. M., Laughlin, G, & Fischer, D. 
2006, ApJ, 637, 1102 

Santos, N. C, Israelian, G., & Mayor, M. 2004, A&A, 415, 1153 

Tereshchenko, V. M. 2001, Astronomy Reports, 45, 1002 

Udry, S., Fischer, D., & Queloz, D. 2007, in Protostars and Planets V ed. 
B. Reipurth, D. Jewitt, & K. Keil, 685-699 

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 

Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 

VandenBerg, D. A., & Clem, J. L. 2003, AJ, 126, 778 

Vogt, S. S. 1987, PASP, 99, 1214 

Vogt, S. S., Allen, S. L., Bigelow, B. C, Bresee, L., Brown, B., Cantrall, 
T., Conrad, A., Couture, M., Delaney, C, Epps, H. W., Hilyard, D., 
Hilyard, D. F., Horn, E., Jern, N., Kanto, D., Keane, M. J., Kibrick, 
R. I., Lewis, J. W., Osborne, J., Pardeilhan, G. H, Pfister, T, Ricketts, 
T., Robinson, L. B., Stover, R. J., Tucker, D., Ward, J., & Wei, M. Z. 1994, 
in Presented at the Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference, Vol. 2198, Proc. SPIE Instrumentation in Astronomy 
VIII, David L. Crawford; Eric R. Craine; Eds., Volume 2198, p. 362, ed. 
D. L. Crawford & E. R. Craine, 362 

Winn, J. N., Henry, G. W., Torres, G., & Holman, M. J. 2008, ApJ, 675, 1531 

Wright, J. T. 2005, PASP, 117, 657 

Wright, J. T, Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261 
Wright, J. T, Marcy, G. W., Butler, R. P., Vogt, S. S., Henry, G. W., Isaacson, 
H, & Howard, A. W. 2008, ApJ, 683, L63 



